Magnetization and Level Statistics at Quantum Hall Liquid-Insulator Transition in 

the Lattice Model 
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Statistics of level spacing and magnetization are studied 
for the phase diagram of the integer quantum Hall effect in a 
2D finite lattice model with Anderson disorder. 
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The way in which the increasing disorder induces the 
insulating state when starting from the integer quantum 
Hall (IQH) state is a topic of controversy between the 
continuum and lattice models of the 2D electronic gas 
in strong magnetic field. The continuum approach pre- 
dicts the crossover between the adjacent quantum Hall 
plateaus, ending up with the insulating state when the 
degree of disorder increases or, equivalently, the limit of 
small magnetic field is considered. This is due to the so- 
called 'floating up' of the critical energies Ec which occurs 
with increasing disorder. {Ec is the energy where the 
localization-delocalization transition takes place in the 
thermodynamic limit). In the critical region, i.e., when 
the Fermi energy Ef crosses an extended state energy Ec, 
the transverse (Hall) conductivity is a!^y — v — 1/2 (at 
the transition between the plateaus v and v — 1)EJ. This 
means that at large disorder (or low field), the cascade 
of transitions must end with v — 1 ^ v = (insulator). 
The experiments give controversial information in what 
concerns the possibility to observe this last transition (see 
Ref.2 and the references therein, Ref.3) The sensible con- 
clusion can be found in Ref.2 and Ref.4 suggesting that 
the theoretical results of the scaling theory, which arc ob- 
tained for zero temperature and infinite systems, cannot 
be checked easily by experiments which are done for finite 
samples and at low (but nevertheless finite) temperature. 
The evaluation of the critical value of the longitudinal 
conductivity is also a difficult task. Lee, Kivelson and 
Zhang show in the frame of corresponding states law that 
^xx — 1/2 for any v; approaching the question in the 
opposite way, Zirnbauer assumes cr^^; = 1/2 and finds 
agreement with the numerical simulations. The numeri- 
cal calculation performed by Huo, Hetzel and Bhatt for 
the lowest Landau level produces also 0.5 ; Huckestein 
and Backhaus show that this value is correct, even in the 
presence of the electron-electron interactiontl. 

More recently, the same problem has been approached 
also in lattice models. The results are again controver- 
sial, since Yang and BhattO, in a discussion based on the 
calculation of the Chern numbers find & tiny floating-up 
of the extended states, while Xie at alQ do not find any 
changing in the position of extended states and affirm 



that the fioating picture is not valid in this model. The 
type of boundary conditions used in the lattice model 
is important; the above mentioned authors use periodic 
conditions. The Dirichlet conditions around the plaque- 
tte, which represent ^a more physical situation, are used 
by Sheng and WengH; their approach, consisting in nu- 
merical calculation of the conductance performed over 
ensembles of disordered finite plaquettes, indicates that 
the QH-insulator transition may occur from any IQH 
plateau of index v directly to the insulator. In the mid- 
dle of the transition region the conductances satisfy the 
relation < Uxx >—< cr^y >— v/2 (here < ... > means 
ensemble average) and the localization length diverges; 
this region is called 'metallic'. 

In this context we study some new relevant features of 
the disordered lattice model in magnetic field with van- 
ishing boundary conditions, the attention being paid es- 
pecially to magnetization, level spacing distribution and 
critical conductances in the metallic region. 

The discussion is based on the spin-less one-electron 
Hamiltonian in perpendicular magnetic field defined on 
a 2D square lattice with N sites in one direction and M 
sites along the other one, which reads as follows: 
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where \n,m) is a set of orthonormal states, localized at 
the sites (n,m), and is the magnetic flux through the 
unit cell measured in quantum flux units. In Eq.(l) the 
hopping integral at = is taken unity serving as the 
energy unit and the diagonal energy £„„ is a random 
variable distributed according to the probability density: 



P(6) 



_{l/W , 
10, 



W/2 < e < W/2 
otherwise, (2) 



(having zero mean e = 0). 

The averaged spectrum of Hamiltonian (1) is depicted 
in Fig.l for (p = 1/10 and the disorder amplitude in the 
range W G [0, 10]. The lines represent the mean eigen- 
values < En > as function of disorder amplitude W. 

In order to study the phase diagram we calculate the 
longitudinal and Hall conductances of this system at con- 
stant magnetic flux and varying disorder. The different 
phases: quantum Hall, metallic, insulating are charac- 
terized not only by conductance but also by the specific 
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distribution of the level spacing and by the current den- 
sity on the plaquette, described by the operator: 

Jnm'' = i tnr^' (r«m -r„',„0 \nm) (n' m' \ +h.c (4) 

(here is the hopping integral between the sites r„m 
and r„/„/). 




FIG. 1. Phase diagram of IQH effect. N=M=10 . The 
insulator regime is depicted for a^x < 0.2 (in /h units) and 
the metal regime for a^x > 0.2 . In the QH regime a^y ~ 
integer and a^x is negligible. The dotted line in the metallic 
region corresponds to the critical points {axx = o^xy)- 

It is opportune to remind previously that for a clean 
system (e„m = 0) with cyclic boundary conditions (i.e 
for a torus) and commensurate values of the magnetic 
flux through the unit cell, the spectrum consists of de- 
generate bands separated by gaps (the well-known Hof- 
stadter butterfly). However, when vanishing boundapj. 
conditions are imposed (i.e., for plaquetteO or cyUndeit2l 
geometry) the gaps get filled with 'edge states', localized 
close to the edges of the sample. The other states, the 
'bulk' ones, remain grouped in bands on the energy scale, 
while geometrically are concentrated in the middle of the 
plaquette. The two types of states differ also by their 
chirality, i.e. by the sign of the derivative dEn/d(j). 
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FIG. 2. The decay of the magnetization (in arbitrary units) 
with increasing disorder for the edge states No. 11 and 10 and 
for the bulk states No. 16 and 12. The distribution of current 
at W = is shown in the insets. <j> — 1/10. 

The effect on the orbital magnetization of each state 
is immediate: the expectation values of the operator 
M = J{r X j(r)) dS calculated on the eigenstates of 



the Hamiltonian (1) have different signs depending on 
whether the state is bulk- or edge-type. 

Fig. 2 shows that the magnetization of the edge eigen- 
state No. 11 is positive Mn > 0, but the bulk eigenstate 
No. 12 has Mi2 < ; the local currents corresponding 
to the two states are also shown in insets. In the same 
figure one anticipates that the increasing disorder pro- 
duces a monotonic decrease of the magnetization. The 
magnetization of all states in the spectrum is shown in 
Fig. 3a and b, for the clean and disordered system, respec- 
tively (the well-known electron-hole symmetry of the Hof- 
stadter spectrum is evident also in the aspect of the mag- 
netization). The disorder effect consisting in the broad- 
ening of the bands and narrowing of the gaps is obvious in 
the second figure. More notable is the magnetization of 
the ground state Mg which can be compared successfully 
with experimental results. Assuming that the spectrum 
is filled up to the Fermi energy Ef, due to the alternat- 
ing sign of Mn in different regions of the spectrum, the 
quantity 

Mg^ J2 (5) 

shows, as function of the number of occupied states, a 
sawtooth aspect (see Fig. 4) which is the same as in the 
experiments by Wiegers at alO, including the fact that 
the jumps of Mg occur at the center of the gaps. In 
our model, the number of teeth depends on the number 
of gaps that can be resolved. For modulated quantum 
wires such a sawtooth aspect |Was obtained also for the 
thermodynamic magnetizatior£3. 




FIG. 3. Magnetization (in arbitrary units) vs. the energ; 
= in fig.Sa and VK = 1 in fig.Sb. N=M=20. (f) = 1/10. 
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no. of electrons 

FIG. 4. The ground state magnetization Mg for N=M=20 
and W = 1 vs. the number of electrons Ne (left scale). Rxx 
and Rxy for the same plaquette coupled to leads (right scale). 
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FIG. 5. Level spacing distribution P{t) for three typi- 
cal situations. Fig.5a gives the distribution of level spac- 
ing for E G [-2.4,-1.3] and W = 3. The dotted line 
is the Wigner-Dyson distribution with B = 2. In fig. 5b 
E e [-2.9, -2.7] and W = 1. The dotted line is the Gaussian 
function whose variance equals the calculated variance of the 
histogram 5{t) = 0.33. In fig.5c E G [-8, -6], W = 13 and 
5{t) = 0.92. The dotted line is the Poisson function. The 
corresponding typical current distributions are shown in the 
insets. N=M=10 and the number of disorder configurations 
is equal to 1000. = 1/10. 



When the Anderson potential (Eq.2) is switched on 
and W is increased continuously, the bands become 
broader and broader, the disorder spreads the states over 
the whole plaquette, giving rise to extended disordered 
states, and finally produces a quasi-continuum of local- 
ized states; even the edge states, which are more ro- 
bust, disappear gradually into the quasi-continuum. The 
nature of the states can be checked by calculating the 
distribution of level spacing for various degrees of disor- 
der, in different domains of the spectrum. Let s„ be the 
level spacing between two consecutive eigenvalues and 
£"„+! and define t„ = Sn/{sn), where (s„) is the mean 
level spacing. 

In Fig. 5 we have three typical distribution functions 
P{t) relevant for different values of disorder and en- 
ergy interval. The Wigner-Dyson (WD) surmise with 
(3 = 2 (unitary) indicates the presence of extended states 
(Fig. 5a) and the Poisson-like distribution (Fig. 5c) shows 
the existence of uncorrelated localized states for strong 
disorder. The case of edge states is also studied, in which 
situation the distribution can be fitted well with the 
Gaussian function (Fig. 5b). Every inset shows a repre- 
sentative distribution of local currents for each case. Very 
illustrative is Fig. 6 which shows the variance 6tn for all 
level spacing of the Hamiltonian (1). One may learn that: 
a) at < 4, for most of the states, 6tn « 0.42, which 
is the typical value for WD distribution with (3 = 2. b) 
for larger W, the variance increases towards 5tn = 1-0 
specific to the Poisson distribution. This value cannot 
practically be reached because of finite dimension of the 
plaquette. c) inbetwcen, at W ~ 6, the variance equals 
0.52, and the probability distribution is well-fitted by the 
WD function with /3 = 1 (so, here the influence of the 
magnetic field is lost and the system behaves 'orthogo- 
nal'), d) the lowest states, originating from the first band 
( n=l,.,6 ) get locahzed faster than the others, while 
the states from the first gap (n=7,.,10) are very robust 
against the localization process. 




'^amplituife of disorder 
FIG. 6. The variance of level spacing 5{tn) (n 



amplitude of disorder W. N=M=10. 



1,..,50) vs. 
The insets show the 



evolution of axx and cjxy at a given number of electrons A'e 
vs. W (for the same plaquette coupled to semiinfinite leads). 
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At last we discuss the transport properties; for this 
purpose the Landauer-Biittiker formaUsm and the tech- 
nics from Ref.9 are used. The conductance in the three 
regimes (IQH, metal and insulator) may be correlated 
with the spectrum characteristics of the isolated system 
discussed above. Since the edge states are responsible 
for IQH effect and they are robust against disorder, this 
regime survives also in the presence of disorder as long 
as W is not too high. In what concerns the transmit- 
tances between different leads, it is well-known that the 
only non- vanishing ones are < Taa+i >j which connect 
consecutive leads (a^lead index) and equal an integer 
value. 

The metallic regime occurs when the transport pro- 
cess takes place on the states which are extended over 
the whole plaquette. In this case also the local current 
is distributed on the whole area allowing for a non-zero 
transmission probability between any pair of leads. This 
behaviour corroborates wit|b-|the multi-fractal properties 
of the local density of statest^. The metallic region can be 
crossed in different ways. Let assume a constant disorder 
(say W =1.0) and change the Fermi level ; in this case, 
metallic regions are intercalated between the QH regions, 
excepting the lowest one which ensures the transition be- 
tween QHl and the insulating phase. The consequence is 
that, when crossing the metallic region, the Hall conduc- 
tance axy will have a smooth decay between consecutive 
plateaus, while axx will differ from zero and have a max- 
imum in the middle of the transition. In Fig. 4, one has 
to notice that Rxx equals Rxy in two places: 1) at the 
transition between QHl and insulator, where one has also 
Cxx « 1/2 in agreement with Huo at alO, and 2) near the 
centre of the spectrum where the system behaves already 
classically. 

Another way to look at the metal-insulator transition 
consists in crossing the metallic zone by increasing the 
amplitude of disorder W, while keeping constant the 
number of electrons (Ne) or the Fermi level. Now the 
transition is of the type — s- and exhibits maxima of 
the longitudinal conductance at a critical disorder Wc 
where, as in Ref.8 , the condition axx = c^xy is fulfilled. 
Such kind of transitions are shown in the insets of Fig. 6 
for Ne=8 (i.e. in the first gap) and for Ne=20 (i.e. in 
the second gap). One may ask what is going on when 
Ne corresponds to a band ; this is the situation for the 
third inset A^e=14. It can be observed that axx is differ- 
ent from zero even at small W, meaning that the bulk 
states are the first ones which become extended under 
the influence of disorder. The dotted line which crosses 
transversally the spectrum in the metallic region of Fig.l 
represents the critical disorder Wc corresponding to each 
N,. 

In the insulating phase, when the Poisson distribution 
of the level spacing is installed, the Hall and longitudi- 
nal conductances tend to zero and the longitudinal re- 
sistance increases exponentially. The transport of the 
electron through the plaquette is performed by tunneling 
on localized states. 



In conclusion, the metallic regime is characterized by 
a Wigner-Dyson distribution of the level spacing with 
f3 — 2 (unitary ensemble). As the system evolves towards 
insulator, the orthogonal WD distribution {(3=1) that 
showes up at a given higher disorder indicates the loss of 
influence of the magnetic field. Simultaneously, the mag- 
netization decays to zero. The QH phase is characterized 
by a Gaussian distribution of level spacing. Due to the 
different chirality of the edge and extended states, when 
crossing the metallic zone (at relatively small disorder) 
the magnetization of the ground state shows a toothsaw 
behavior as function of the filling factor similar to exper- 
imental results. 
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